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Abstract 

Using higher order intertwining operators we obtain new exactly solvable potentials admitting quasi- 
normal mode (QNM) solutions of the Klein-Gordon equation. It is also shown that different potentials 
exhibiting QNM's can be related through nonlinear supersymmetry. 
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1 Introduction 



Quasinormal modes (QNM) are basically discrete complex freequency solutions of real potentials. They 
appear in the study of black holes and in recent years they have been widely studied pQ . Interestingly 
QNM's have also been found in nonrelativistic systems [2J. However as in the case of bound states 
or normal modes (NM) there are not many exactly solvable potentials admitting QNM solutions and 
often QNM frequencies have to be determined numerically or other approximating techniques such as the 
WKB method, phase integral method etc. Consequently it is of interest to obtain new exactly solvable 
potentials admitting such solutions. 

In the case of NM's or scattering problems a number of methods based on intertwining technique 
e.g, Darboux algorithm [3J, supersymmetric quantum mechanics (SUSYQM) [1] etc. have been used 
successfully to construct new solvable potentials. Usually the intertwining operators are constructed using 
first order differential operators. However in recent years intertwining operators have been generalised 
to higher orders [3 [9] and this has opened up new possibilities to construct a whole new class of 
potentials having nonlinear symmetry. In particular use of higher order intertwining operator or higher 
order Darboux algorithm leads to nonlinear supersymmetry. 

QNM's are associated with outgoing wave like behaviour at spatial infinity and unlike normal modes 
(NM's), the QNM wave functions have rather unusual characteristic (for example, wave functions diverg- 
ing at both or one infinity) [10. . Such open systems have been studied using (first order) intertwining 
technique [11] . Recently it has also been shown that open systems can be described within the framework 
of first order SUSY [10j . Here our objective is to examine whether or not intertwining method based on 
higher order differential operators can be applied to open systems. For the sake of simplicity we shall 
confine ourselves to second order intertwining operators (second order Darboux formalism) and it will 
be shown that the second order Darboux algorithm can indeed be applied to models admitting QNM's 
although not exactly in the same way as in the case of NM's. In particular we shall use the second order 
intertwining operator to the Poschl- Teller potential to construct several new potentials admitting QNM 
solutions. It will also be shown that such potentials may be related to the Poschl- Teller potential by 
second order SUSY. The organisation of the paper is as follows: in section 2, we present construction of 
new potentials using second order intertwining operators; in section 3, nonlinear SUSY underlying the 
potentials is shown and finally section 4 is devoted to a conclusion. 



2 Second order intertwining approach to quasinormal modes 

Two Hamiltonians H and Hi is said to be intertwined by an operator L if 

LH = HiL (1) 

Clearly if ip is an eigenfunction Ho with eigenvalue E then Lip is an eigenfunction of Hi with the same 
eigenvalue provided Lip satisfies required boundary conditions. If L is constructed using first order 
differential operators then intertwining method is equivalent to Darboux formalism or SUSYQM. In 

particular if Vo is the starting potential and L = — — \-W(x), then the isospectral potential is V\ — 
dW 

Vq + 2— — HI. Similarly if L is generalised to higher orders then it is equivalent to higher order Darboux 
ax 

algorithm or higher order SUSY. 
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Let us now consider L to be a second order differential operator of the form [7] 
L 
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where ipi and tpj are eigenfunctions of Hq corresponding to the eigenvalues wf and us 2 and Wij = 
(ipiip'j ~ ifypj) is the corresponding Wronskian. Then the isospectral partner potential V^x) obtained via 
second order Darboux formalsim is given by 

V 2 (x) = V (x) 



dx'- 



■ log Wi tj (x) 



The wave functions ipi( x ) an d </>i(^) corresponding to Vo(x) and V2{x) are connected by 



<j> k {x) = Lip k (x) 
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The eigenfunctions obtained from ip^ and ipj are given by 



f(x) oc 



W id (x) 



g(x) oc 
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It may be noted that in the case of normal modes, the new potential would be free of any new singularities 
if the Wronskian Wij (x) is nodeless. This in turn requires that the Wronskian be constructed with the 
help of consecutive eigenfunctions (i.e, j = i + 1). Also, the eigenfunctions f(x),g(x) in ([5]) are not 
acceptable because they do not satisfy the boundary conditions for the normal modes and in any case 
they are not SUSY partners of the corresponding states in the original potential. Thus in the case of 
normal modes the spectrum of the new potential is exactly the same as the starting potential except for 
the levels used in the construction of the Wronskian. However we shall find later that not all of these 
results always hold in the case of QNM's. 



Let us now consider one dimensional Klein- Gordon equation of the form [lOj 

{d 2 ~d 2 x + V(xM(x 1 t)=Q 



The corresponding eigenvalue equation reads 

Hip n = <jj%ip n (x), H 



dx 2 



V{x) 



(6) 



(7) 



The QNM solutions of the equation {7J are characterised by the fact that they are either (1) increasing 
at both ends (II) (2) increasing at one end and decreasing at the other (ID,DI). The wave functions 
decreasing at both ends (DD) correspond to bound states or NM's. In the case of QNM's the eigenvalues 
(u>n) m ay be complex or real and negative. If Re(uj n ) ^ then the SUSY formalsim can not be applied 
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since in that case the superpotential W(x) becomes complex and consequently one of the partner potential 
becomes complex. So we shall confine ourselves to the case when Re(u) n ) = i.e, lo 2 ^ are real and negative. 
We would also like to mention that in case Eq.([7]) is to be interpreted as a Schrodinger equation one just 
has to consider the replacement io 2 t — > E n . 

There are a number of potentials which exhibit QNM's. A potential in this category is the inverted 
Poschl- Teller potential. This potential is used as a good approximation in the study of Schwarzschild 
black hole and it is given by 

Vq{x) = vsech 2 x (8) 

Eq.([7]) for the potential © can be solved in different ways. One of the simplest way is to apply the shape 
invariance criteria [4] and the solutions are found to be [12j [T3l [14] 

w± = -i(n-A±), A± = - l -±^\~v = - l -± q (9) 

tf(x) = (sechx)^-^^ + q \-q- 1 - iw± 1 + * nhx ) (1Q) 

We note that the behaviour of the wave functions (TTTJ)) i.e, whether they represent a NM or QNM depends 
on the value of the parameter v. For v S (0,1/4) i.e, A + G (— 1/2,0), A~ G (—1, —1/2), the wave functions 
represent outgoing waves and are of the type (II). For v < (i.e, A + > 0, A~ < 0) then the wave functions 
represent NM's when n < A + while for n > A + they are QNM's. On the other hand the wave functions 
are always QNM's corresponding to u>~. It may be noted that for the QNM's the wave functions (TTU)) 
for even n are nodeless while those for odd n have exactly one node at the orign. This behaviour of the 
wave functions is quite different from those occuring in the case of NM's. Using the procedure mentioned 
above we shall now construct new exactly solvable potentials admitting QNM solutions. 



2.1 Construction of isospectral partner potential using NM's 

Case 1. Vq with two NM's: In order to apply the second order intertwining approach one may start 
with a potential Vq{x) admitting (1) at least two NM's and the rest QNM's or (2) only QNM's. We begin 
with the first possibility. Thus we consider v — —5.04 so that there are two NM's. In this case we obtain 
from ([9|) A ± = 1.8, -2.8. Thus the NM's correspond to lu£ = 1.8i,Lof = O.Si and are given by 

^+(x) = (sech x) A+ , iPt{x) = (sech x)^ A+ - 1 ^ 2 F 1 (2A + , -1, A+, 1 + ta ^ hx ) (n) 

The QNM's in this case correspond to the frequencies w^,n = 2,3,... and are given by ip^(x). Also 
for lu~ , n — 0, 1, 2, ... there is another set of QNM's and the corresponding wave functions are given by 
tp~(x). In this case A~ < and consequently there is no NM. 

Let us construct a potential isospectral to ([5]) using the NM frequencies Uq and u)± . Then from (TlT))) 
the Wronskian Wq 1 is found to be 

W + = -(sechx) 2 ^- 1 (12) 

Clearly W^~ x does not have a zero. In this case the new potential V^~(x) is free of singularities and is 
given by 

d 2 

V+(x) = V (x) - 2— logWo+^x) = -(1 - A+)(2 - A+)sech 2 x (13) 
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Using the value of A + it is easy to see that the new potential V 2 (x) in (fT5|) does not support any bound 
state but only QNM's. This is also reflected by the explicit expressions for the wave functions. Using 
we find 

f+(x) = (sechx) {1 - A+) , g + {x) = ~{sechx)- A+ tanhx (14) 

From (fT4"|) it follows that the above wave functions are QNM's corresponding to — ui^ and — uJq respec- 
tively. Note that these two QNM's are new and were not present in the original potential. This in fact 
is where the behaviour of the new potential is different from the usual case. In the case of potentials 
supporting only NM's the wave functions f + (x),g + (x) obtained through §5§ do not have acceptable be- 
haviour. However in the present case both these wave functions become QNM's instead of NM's and 
they have acceptable behaviour at ±00 as can be seen from (1141) as well as from figure 1. The other wave 
functions <^>+(x),n = 2,3, ... corresponding to QNM frequencies w+ = —i(n — A + ) can be obtained using 
(j4j and are given by 

4>+(x) = (sechx)( A [(n - l)n F n tanh 2 x 



+ci(2n — 3) F n+ i sech 2 x tanhx 



(15) 



where 



ci 



+ cic 2 F n+2 sech 4 x] , n = 2, 3, 



n(2A+ -n+ 1) 
' 2{A+ -n + 1) 



C2 = 



{-n+l)(2A+ -n + 2) 
2{A+ - n + 2) 



F„ 

F„+i 



iFi(-n,2A+ -n + 1,A+ -7i+ 1, 



l-\-tanhx 



iF x (-71 + 1, 2A+ - n + 2, A+ - n + 2, l+lfih2L) 



(16) 



F n+2 = 2 Fi{-n + 2,2A+ -n + 3,A+ 



\+tanhx 
2 



To see the nature of the wave functions (TT5|) we have plotted 0^~(x) and (j>t{x) in figure 1. From the 
figure it can be seen that these wave functions are indeed QNM's and for even n they do not have nodes 
while for odd n they have one node at the origin. We would like to point out that the new potential 
V 2 + (x) has two more QNM's than Vo(x). Thus except for two additional QNM's, the QNM frequencies 
w+ is common to both Vq(x) and V 2 + (x). We now examine the second set of solutions corresponding to 



It can be shown by direct calculation that the new potential (jT3|) also possess this set of solutions. 

Case 2. Vq(x) with three NM's: Let us now consider the potential ([5]) supporting three NM's. A 
convenient choice of the parameter is v = —6.2 so that A + = 2.04, A~ = —3.04. We shall now construct 
the new potential using the NM frequencies uj^ = 1.04?, and lo 2 = 0.042.. The Wronskian W^ 2 is found 
to be 



W+(x) = 



(sechx) 



2A^ 



2(A+ - 1) 



■[A+-2- {A + - l)cosh2x] 



(17) 



Now using © we obtain 



V 2 + (x) = -(A+ - 1)(A+ - 2)sech 2 x + 8(A + - 1) 



(A+ - 2)cosh2x - {A+ - 1) 
p+ - l)cosh2x - (A+ - 2)p 



(18) 



To get an idea of the potential, we have plotted V 2 + (x) in figure 2. From figure 2, it is clear that V 2 + (x) 
supports at least one NM. Next to examine the wave functions we first consider f + (x) and g + {x). ^From 
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the relation © we obtain 

+ _ 2(A+ - l)(sechx)- A+ tanhx + _ (sechx)-( A++1 ^ [1 - (2A+ - l)tanh 2 x] 
* (X) ~ (A+ - l)cosh2x - (A+ - 2) ' ff ^ ~ (1+ - l)cos/i2a; - (A+ - 2) (19) 

Also from ((4]) it follows that 

00 W_ (A+-l)c OS / l 2a; + A+-2 ( j 

^From (fl9|) it follows that f + (x) and g + (x) are new QNM's corresponding to frequencies — uf^ = — 0.04i 
and — uj^ = — 1.04z respectively. The former has one node the later has two nodes. The nodal structure 
of the QNM wave functions are different from those obtained earlier. The reason for this is that since we 
started with the first and second excited state NM's and the Wronskian Wj 2 is nodeless, the behaviour 
of the original wave functions tpi 2 { x ) are retained by f + (x) and g + {x). However, 4>q{x) is a NM at 
u)q = 2.04:2 and it does not have a node because i/jq(x) does not have one. Also other QNM wave 
functions cj>+(x),n = 3,4, ... have either no node or one node. In figure 3 we have plotted the some of 
the wave functions. We also note that although the potential in (p~3j) is of a similar nature as J5]), the 
potential (I18p is of a completely different type. In particular it is a non shape invariant potential. Finally 
we discuss the possibility of a second set of solutions for the potential (fT8|) . We recall that the existence 
of two sets of solutions for the potential JHJ (or (fl3|) ) was due to the fact that the parameter v could be 
expressed as a product of two different parameters A^. However in the case of (TT5|) the entire potential 
can not be expressed in terms of two distinct parameters because of the presence of the second term. 
Consequently the potential (fT8|) has only one set of solution mentioned above. 



2.2 Construction of isospectral partner potential using QNM's 

Case 1. Potential based on consecutive QNM's: Here we shall construct isospectral partner of a 
potential which has only QNM's. Thus we consider v = 0.24 and in this case A^ = —0.4, —0.6. We 
consider the A + sector and begin with the freequencies Wq and uif . In this case the expression for 
the Wronskian W^ 1: the new potential V^{x) and the QNM wave functions can be derived from the 
expressions (|12[) . (|l3j) and (fl"5|) respectively except that we now have to use a different parameter value. 
Thus the new potential is given by 

V^(x) = -3.36 sech 2 x (21) 
For this potential the NM's corresponding to — uiq = OAi and —u>x = 1.4i are given respectively by 

/+ = (sechx) 1A , g + = (sechxf Hanhx (22) 

Clearly these NM's are not SUSY partner of any levels in H . The QNM's are correspond to w+ = 
— i(n + 0.4), n = 2,3, ... and are given by (fTS)) with A + = —0.4. We have plotted some of the wave 
functions in fig 4. From the figure we find that the wave functions f + (x) and g + (x) correspond to NM's 
and the other wave functions represent QNM's which are the SUSY partners of the QNM's in Hq. We 
note that as in (fT3"l) the potential (f2Tj) has two sets of QNM's, the second of which corresponds to u>~. 



Case 2. Potential based on non consecutive QNM's: Here we shall consider the previous pa- 
rameter values (i.e, A + = —0.4) and construct the new potential using the non consecutive levels Uq and 
ui^. In this case the Wronskian is given by 



( sechx)^ ~ 3 ) 

Ks = 2(A+-2) [(9 ~ 6A )Unh X + 3] (23) 
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It can be shown that the Wronskian (f2"3")) is nodeless. Now using the © the new potential is found to be 

(A+ - 2)[2{A+(A+ - 2){3A+ - 7) - 2A+(A+ - 1)(A+ - 4) cosh2x - (3 - 2A+) 2 (A+ - 1) sech 2 x] 



V 2 + (x) 



[1-A+ + (A+ - 2) 

(24) 

The potential (|24|) is free of any singularity and is plotted in fig 5. From figure 5 we find that it supports 
NM's. Also as explained earlier, this potential has also one set of solution. We now consider the wave 
functions corresponding to i/jq(x) and ip^(x). These are obtained from §5§ and are given by 

f+( x ) - 2(A+-2) ( sechx )(3-A + ) a+(x) _ (1 - 2A+)tanh 2 x + 3 ( sechx )d-A + ) 

J W " (9 _ 6A+)tanh 2 x + 3 ( sechx > > 9 W - (9 _ 6A +) tanh 2 x + 3 smhx ( secllx > 

(25) 

The above wave functions (with zero and one node respectively) represent NM's corresponding to — = 
3Ai and — Wq = OAi. The other wave functions can be obtained through ((4]). The two QNM wave 
functions lying between ljq and tu^ are 4'i2(. x ) corresponding to cu^ = — lAi and u>2 = — 2Ai. We have 
plotted these wave functions in fig 6. From figure 6, it can be seen that f + {x) and g + (x) are NM's while 
(jy[ 2 (x) are QNM's, with the later having two nodes. The rest of the QNM wave functions corresponding 
to the freequencies oj+ = — (n + 0.4)i, 0,3 are given by <ft£(x) and they have either zero or one node. 



3 Polynomial SUSY 

In first order SUSY, the anticommutator {Q, Q'} of the supercharges is a linear function of the Hamil- 
tonian. On the other hand in higher order SUSY, {Q, Q^} is a nonlinear function of the Hamiltonian. It 
will be shown here that the Hamiltonians Hq and Hi are related by second order SUSY. To this end we 
define the supercharges Q and as follows: 



^ L / ' V 

where the operator L is given by ([2]). 

Clearly the supercharges Q and are nilpotent. We now define a super Hamiltonian H of the form 

I) < 27 ' 

It can be easily verified that Q, and H satisfy the following relations : 

[Q,H] = [Q\H]=0 (28) 
Then the anticommutator of the supercharges Q and Q* is given by a second order polynomial in H : 

L^L \ /„ <5 x2 



H S s = {Q,Qi}=^ Q L " Lt j = [ H +2) - cI < 29 ) 
where X is the 2x2 unit matrix and 

6 = -{tf+ufi , c =(^-^l) (30) 



7 



Also we have 

[Q,H„] = [Q\H a8 ]=0 (31) 

The relations (f2"9"| and (j3"Tj) constitute second order SUSY algebra. 

As an example let us consider the potentials ||5J) and (fT3"|) . The corresponding Hamiltonians f/o and 
H2 are obtained from ([7]). In this case 5 = 0.5416 and c = 0.2916 so that from (|29|) we get 

H ss = (i? + 0.5416) 2 - 0.29162: (32) 

In a similar fashion one may obtain H ss for the other pair of potentials. 



4 Conclusion 

Here we applied the second order Darboux algorithm to the Poschl- Teller potential and obtained new 
exactly solvable potentials admitting QNM solutions. We have considered a number of possibilties to 
construct the new potentials e.g, starting from NM's or starting from QNM's. It has also been shown 
that the new potentials are related to the original one by second order SUSY. We feel it would also be 
also useful to analyse the construction of potentials using various levels as well as for different values of 
the parameter v (for example, v = half- integer) [lOj . Finally we belive it would be interesting to extend 
the present approach to other effective potentials appearing in the study of Reissner-Nordstrom, Kerr 
black hole etc. 
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Fig 1. Graph of f(x), g + (x), (J) 2 + (x) and (|)3 + (x) 
for A + = 1.8 
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Fig 2 . Graph of V 2 + (x) for A + = 2.04 
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Fig 3. Graph of f^x), g + (x), (j) + (x) and <j)3 + (x) 
for A + = 2.04 
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Fig 4. Graph of f(x), g + (x), ())2 + (x) and ()>3 + (x) 
forA + = -0.4 
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Fig 5 . Graph of V 2 + (x) for A + = - 0.4 
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Fig 6. Graph of f^x), g + (x), (j)i + (x) and (j)2 + (x) 
for A + = - 0.4 



